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Many-body calculations of the total energy of interacting Dirac electrons in finite graphene sam- 
ples exhibit joint occurrence of cusps at angular momenta corresponding to fractional fillings charac- 
teristic of formation of incompressible (gapped) correlated states (y = 1/3 in particular) and opening 
of an insulating energy gap (that increases with the magnetic field) at the Dirac point, in corre- 
spondence with experiments. Single-particle basis functions obeying the zigzag boundary condition 
at the sample edge are employed in exact diagonalization of the interelectron Coulomb interaction, 
showing, at all sizes, mixed equal-weight bulk and edge components. The consequent depletion of 
the bulk electron density attenuates the fractional-quantum-Hall-effect excitation energies and the 
edge charge accumulation results in a gap in the many-body spectrum. 
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I. INTRODUCTION 

The isolation of graphene sheets^ was soon followed by 
experiments^ on the anomalous integer quantum Hall 
effect (IQHE), demonstrating the formation of Landau 
levels, described with non-interacting electron theory. 
Subsequently, theoretical explorations 4,5 appeared per- 
taining to the fractional quantum Hall effect (FQHE) in 
graphene, which, in contrast to the IQHE, is known to 
be the manifestation of collective correlated interacting 
electrons phenomena. For Dirac electrons in the low- 
est (n = 0) Landau level (LLL), it was found^' 5 (in a 
spherical geometry) that the behavior of the incompress- 
ible states in boundless graphene models emulates that 
of the well-known FQHE in the semiconductor 2D elec- 
tron gas£~— Recent experiments^^ employing small (1 
fj,m) suspended single-sheet samples have shown that the 
appearance of the FQHE in graphene is accompanied by 
the emergence of an insulating phase. The physics that 
underlies this joint occurrence, coupled with the system- 
atic theoretical overestimation 4 ' 5 of the measured FQHE 
excitation energies, remains unresolved; similar findings 
in recent experiments on bilayer graphene^i further high- 
light these open issues. 

In this paper, we show, using exact diagonalization 7 
(EXD) (in the quantum Hall regime) of the many-body 
hamiltonian for interacting Dirac electrons occupying the 
LLL of finite planar graphene samples, that the incom- 
pressible correlated electron physics is modified [using 
zigzag boundary conditions 12 (ZBC)] in two main ways: 
(i) the FQHE excitation energies are significantly atten- 
uated due to depletion of the density of the bulk compo- 
nent of the Dirac electrons spinor, and (ii) the concurrent 
accumulation of charge at the sample edge underlies the 
opening of a gap in the many-electron spectrum (increas- 
ing with the magnetic field, as observed experimentally), 
relating to the aforementioned insulating phase. 

Underlying the above many-body behavior is the ex- 



istence of LLL single-particle states of mixed charac- 
ter, exhibiting bulk- [Darwin-Fock-type 7 (DF)] and edge- 
type components of equal weights^ (independent of sam- 
ple size), each residing on a different graphene sub- 
lattice. With interelectron repulsion, the bulk com- 
ponent exhibits FQHE characteristics (with excitation 
energy gaps reduced by a factor of four compared to 
boundless graphene). The edge component is associ- 
ated with formation of a single-ring rotating Wigner 
molecule 7 (RWM). The size-independent influence of the 
graphene edge on the many-body properties ushers a new 
paradigm, contrasting the commonly accepted expecta- 
tion that materials bulk properties are unaffected by the 
detailed configuration or conditions at the surface, 14 This 
unique behavior derives from the two-sublattice topol- 
ogy of the graphene net modeled here by the continuous 
Dirac- Weyl (DW) equation* 1 ^ 

An important clement in our approach is formulation 
of the solutions of the Dirac- Weyl equation (with the 
zigzag boundary conditions) using the Kummer conflu- 
ent hypergeometric function. 113 This formulation provides 
systematic analytical insights into the nature of the DW 
single-particle states (bulk and edge), as well as it per- 
mits efficient and accurate numerical evaluations (with 
the help of algebraic computer languages^). 

In Sec. we describe the confluent-hypergeometric- 
function formulation of the solutions of the DW equation 
and construct the LLL DW spinor in terms of equal- 
weight bulk and edge components. Details of the solu- 
tions of the DW equation in polar coordinates are given 
in Appendix A. 

A brief description of the exact-diagonalization many- 
body method, and the contributions of the bulk and edge 
spinor components to the Coulomb interaction matrix 
elements are given in Sec. Mil 

In Sec. IIVI we display results for the ground-state and 
excited spectra, and their relation to the v = 1/3 FQHE 
incompressible state and to the insulating behavior at 
v = 0. 
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A summary is presented in Sec. El 

Finally, in Appendix B we distinguish the mixed bulk- 
edge Dirac spinor states (forming the quasidegenerate 
LLL manifold in finite graphene samples with ZBC), 
where one of the spinor components is located in the bulk 
of the sample and the other localized at the edge, from 
the "double-edge" states where both components of the 
spinor represent edge states. In the "double-edge" case, 
one component corresponds to a bulk orbital which trans- 
forms (when the orbital's centroid falls near the graphene 
sample boundary) into an edge state familiar from the 
theory-- of the nonrelativistic integer QHE. The other 
component of the "double-edge" spinor, as well as the 
edge orbital of the aforementioned mixed bulk-edge one, 
corresponds to an edge state unique to the two sublat- 
tice topology of graphene. We discuss and illustrate the 
natural way in which the above classification of Dirac- 
electron edge states arises within the framework of the 
Kummer function formalism for the solution of the DW 
equation described in Sec. Inland Appendix A. 



II. THE SINGLE-PARTICLE 
LOWEST-LANDAU-LEVEL MANIFOLD WITH 
ZIGZAG BOUNDARY CONDITIONS 

We first discuss the solution of the Dirac- Weyl equa- 
tion in polar coordinates under the imposition of the 
zigzag boundary condition. We model the low-energy 
noninteracting graphene electrons (around a given K 
point) via the continuous DW equation^ Circular sym- 
metry leads to conservation of the total pseudospin 12 
j = I + a z , where / is the angular momentum of a 
Dirac electron. As a result, we seek solutions for the 
two components ip A (r) and 4> B ( r ) (associated with the 
two graphene sublattices A and B) of the single-particle 
electron orbital (a spinor) that have the following general 
form in polar coordinates: 
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The angular momentum takes integer values; for sim- 
plicity in Eq. ^ and in the following, the subscript I is 
omitted in the sublattice components tp A , ip B and \ A L > 
X B 

With Eq. ([1]) and a constant magnetic field B (symmet- 
ric gauge), the DW equation reduces (for the K valley) 
to 



Tx x ~~ x [ l + Y^' 



(2) 



where the reduced radial coordinate x — r/ls with 
Ib = y/tic/ (eB) the magnetic length. The reduced 
single-particle eigenenergies e = E/(Tivf/Ib), with vf 
the Fermi velocity. Since the properties of the solutions 



of the DW equation with a magnetic field using the ZBC 
are not widely known, we outline pertinent details on this 
subject in the following (see also Appendices A and B). 

For a finite circular graphene sample of radius R, we 
seek solutions of Eq. @ for e ^ with I < — 1 (cor- 
responding to j < —1/2) that are regular at the origin 
(x = 0). The general form of solution is 
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where M(a, 6, z) is Kummer's confluent hypergeometric 
function^ and C is a normalization constant. 

The ZBC requires that one sublattice component van- 
ishes at the physical edge situated at the finite radius R 
of the graphene sample. Therefore at least one of the 
Kummer functions in Eq. ([3]) must exhibit zeros (nodes) 
for x > 0. The search for such zeros is greatly facilitated 
with the use of the following theorem*^ For the Kum- 
mer function M(a, b, z) (z being real), it is known that, 
if b > [which is the case for both Kummer functions 
in Eq. (J3J], there are no zeros of M(a,b,z) for z > 
if a > 0, and there are precisely — \_a\ zeros if a < 0; 
the floor function is defined as: [ a J = ?i if and only if 
n < a < n + 1, where n is any positive or negative inte- 
ger (including zero). 

The \ B LLL component does exhibit precisely one 
zero, enabling a transcendental equation (for a given xr) 
for the single-particle energies e 



X B (e,l,x R ) =0, or M 



£ 

~~2 



= 0. 



(4) 



The solutions to Eq. (j4j (for a range of xr » 1, i.e., 
high magnetic fields, and for any I < — 1) give states 
(which compose the LLL manifold) with near vanish- 
ing energies < e/%/2 << 1. Consequently, within the 
graphene sample, when R/Ib » 1, X B can be approxi- 
mated by a bulk- type DF orbital, i.e., 



\b) =Ce-^l\\ l \~\ 



(5) 



with C = 1/ ^/2l / l _1 r(|Z|) normalizing \b) to unity in 
[0,xr]. 

From the theorem discussed above, it follows that the 
X A LLL component cannot have any zeros, and thus for 
Xr = R/Ib >> 1 it will develop into an edge state 

\e) = -|C|e* 2 /VI, (6) 

due to the asymptotic behavio r 15 ! 19 

M(l - e 2 /2, |*| + 1, x 2 /2) - exp(a; 2 /2)a; fl 21 ' 1 (7) 

for x xr. The normalization constant C conforms to 
the requirement that x A and x B [ see Eq. (U}] should have 
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FIG. 1: EXD energies for N = 6 Dirac electrons. Ener- 
gies for the lowest (yrast, solid dots) and first-excited states 
(solid squares) are plotted as a function of the total angular 

momentum M = X^i=i mi ( mi = !'*!)• "b"> " e '\ an d "be" 
denote, respectively, the bulk, edge, and cross bulk-edge con- 
tributions to the total energies. The "b" contribution reflects 
a 50% depletion of the bulk component due to the form of 
the Dirac-electron spinor in Eq. (|8]). The arrow in the "b" 
curve points to the v — 1/3 (Af = 45) excitation energy gap 
(the difference between first-excited and yrast energies). The 
radius of the sample was taken as R = 301b- k = 1 for sus- 
pended graphene. 
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FIG. 2: EXD energies for N = 7 Dirac electrons. Energies 
for the lowest (yrast, solid dots) and first-excited states (solid 
squares) are plotted as a function of the total angular momen- 
tum M = Yli=i mi i m i = IM)- "b", " e ", an d "be" denote, 
respectively, the bulk, edge, and cross bulk-edge contribu- 
tions to the total energies. The "b" contribution reflects a 
50% depletion of the bulk component due to the form of the 
Dirac-electron spinor in Eq. J5J. The radius of the sample 
was taken as R — 251b- k = 1 for suspended graphene. 



III. EXACT DIAGONALIZATION OF THE 
MANY-BODY PROBLEM 



equal weights for the e's obtained via solution of Eq. 
C is determined from C 2 J^ R e x / 2 x 2 ^ +1 dx = 1. 

Finally, for R/Ib >> 1, the following simplified ex- 
pression for the LLL spinor can be written: 



i 

V2 



e l ^|e)/V27T 
s <(i+i)*|6)/^ 



(8) 



In the EXD method^ the many-body wave function is 
given as a linear superposition 

mi + ...+mjv — M 

V(m 1> ...,m N )^-' mN [r 1 ,...,r N ], (9) 

mi<...<mjv 

where the Slater determinants \j/ m ^---. m Jv are constructed 
out of the quasidegenerate (e ~ 0) LLL spinor orbitals 
ipl(r) [see Eq. (|5j)], where we define m = \l\; I < — 1. 
We omit the Coulomb coupling between valleys, and also 
consider fully spin-polarized electrons. 

The eigenvectors V's in Eq. © and corresponding 
many-body energies _E EXD are calculated through a di- 
rect matrix diagonalization of the many-body hamilto- 



where \b) and |e) are each normalized to unity in [0,£r]. 
A numerical graphical illustration showing that the gen- 
eral expression for the spinor components in Eq. Q can 
be well approximated by Eq. (|5J| is given in Appendix B. 



H = J2v c (i,j), 



i<3 



(10) 



where Vc(i,j) = e 2 /(K|rj — i\,|) is the two-body 
Coulomb interaction, within the Hilbert space of Slater 
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determinants^ <J/ mi ' •••> m «_ B e i n g constant (here ss 0), the 
kinetic energy terms in the many-body hamiltonian Eq. 
([TOf have been neglected.— ^ 

The EXD computations require an accurate evaluation 
of the two-body Coulomb matrix elements 



dridr 2 ip* mi {v x )V4 2 ( r 2 ) V c (1 , 2)ip„ l3 (r x )ip mi (r 2 ) , 

(11) 

which expand to a sum of four similar integrals. Denoting 
\b m ) = ie^ +1 ^\b)/y^ and \e m ) = e fl *|e)/v^7r, this 
expansion is written as 

^((foi&2|^3&4) + (eie 2 |e 3 e 4 ) + (6ie 2 |& 3 e 4 ) + (ei& 2 |e 3 6 4 )). 

(12) 

Note that, due to the equal weights of the bulk-like and 
edge-like components, a prefactor of 1/4 appears in front 
of each term in Eq. (IT21) ; the consequences of this pref- 
actor are discussed in Sec. IIV Al 



IV. RESULTS 

A. Energetics 

In Figs. [T] and [U we display total EXD energies for the 
yrast (i.e., the lowest in energy) and first-excited states of 
spin-polarized electrons in a graphene finite sample with 
N = 6 and N — 7 Dirac electrons as a function of the 
total angular momentum M (which is conserved in the 
EXD calculation). Moreover, the displayed energies cor- 
respond to: (i) only the bulk- type terms ("b"), (ii) com- 
bined bulk- type and edge-type terms ("b+e"), and (iii) 
with the inclusion of cross bulk-edge terms ( "b+e+be" ) 
in the two-body Coulomb matrix elements [see Eq. (fT2|)]. 
Each of the terms ("b", "e", and "be") is reduced by a 
factor of 1/4 as discussed above. For all cases, we note the 
appearance of cusp states (states with enhanced stabil- 
ity and larger excitation gaps relative to their immediate 
neig hborhoodZ&2fi) at the magic angular momenta 

M m = M + kN or M m = M + k(N - 1), (13) 

where M = N(N ~ l)/2, k = 0,1,2,3,..., and N is 
the number of electrons. For the bulk component, the 
cusp states are interpreted^^ as a signature of the forma- 
tion of correlated incompressible states that may exhibit 
certain correlated liquid characteristics and underlie the 
FQHE physics in the thermodynamic limit. In particu- 
lar, the magic angular momenta are associated^^ with 
fractional fillings v — Mo/M m . Note that the angular 
momenta of the Laughlin FQHE function^ correspond to 
Eq. (p~3|) with k = (N — l)p for the first case and k = Np 
for the second case, with the associated fractions given 
by v = l/(2p + 1) (p being a positive integer). 

The EXD energies displayed in Fig. [T] and [5] [and ad- 
ditional calculated results for other sizes N = 3 — 5 and 
8 (not shown)] reveal three prominent trends: 
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FIG. 3: Variation of the insulating gap fx. (a) As a function 
of B when R = 1 /im. (b) As a function of R when 5 = 8 
T. Note the fi oc B 054 increase shown in (a) and the effective 
saturation for large R exhibited in (b). a = 2. 



(I) The "b+e" energies are always higher than those 
calculated only with the bulk-type component (denoted 
as "b"). More importantly (as checked for various sam- 
ple sizes, e.g., R/Ib = 15,25 and 30), the difference be- 
tween the "b+e" (i.e., including the edge contribution) 
and "b" EXD energies (for any N) approaches asymp- 
totically for M -> oo the value AE t (N)/4; AE^N) is 
the electrostatic energy of N point charges localized on 
the perimeter at the vertices of a regular polygon, i.e., 



AEx(N) 



NS N e 2 
4xr kU 



(14) 



with S N = J2f= 2 ( S HU - iW^I) L - Th e factor 1/4 
is due to the half weight of the electrons accumulating 
on the physical edge [recall the spinor in Eq. ((5J]. The 
result in Eq. (fT4| (deduced from analysis of the EXD 
computed energies, see, e.g., Figs. Q] and [2]) suggests that 
a single-ring Wigner molecule is formed at the edge; this 
is in agreement with the CPD analysis in Sec. IIVBI 

(II) The "b+e+be" energies are always higher than 
the "b+e" ones. Moreover, the difference between the 
"b+e+be" and "b+e" EXD energies (for any N) ap- 
proach asymptotically for M — > oo the value AE 2 (N) /4, 
where Ai? 2 (iV) is the electrostatic energy due to the re- 
pulsion between a point charge of strength Ne located at 
the center and A^e charges distributed at the perimeter 
of the sample, 



AE 2 (N) 



N(N - 1) e 2 



IR 



kIj 



(15) 



The factor 1 /4 is again due to half of the electrons being 
accumulated on the physical edge, while the other half 
being located in the bulk component; the factor A^ — 1 
accounts for the self-interaction correction. 

It is apparent that the EXD density of states associated 
with the "b+e+be" LLL spectrum will show (for a given 
N) the opening of an interaction-induced energy gap 



A£(A0 = AE 1 (N)/A + AE 2 {N)/A 



(16) 



compared to that associated with the bulk-only ("b") 
spectrum. The emergence of this edge-induced, electro- 
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FIG. 4: CPD for the edge component portraying the forma- 
tion of a single-ring RWM. An illustrative yrast cusp state 
(with N = 6 and M = 39) in a graphene sample of radius 
R/Ib = 30 is considered. The azimuthal CPD at the phys- 
ical edge of the sample is plotted, i.e., V(r, ro) as a func- 
tion of r = 30(cos 6, sin 6). The fixed point is located at 
ro = 30(cos 77, sin n). Note the five maxima corresponding 
to the remaining five electrons, with the sixth one taken at 
ro . r and ro are in units of Is ■ 



static energy gap is related to (see below) the experi- 
mentally observed2J£ insulating behavior at the charge 
neutrality point {v = 0). 

(Ill) Due to the fact that both the "e" and "e+be" 
contributions (expressed in units of e 2 /kIb) decrease as 
a function of R/Ib [see points (I) and (II) above, and 
the factor l/x R = Ib/R in Eqs. ([14]) and (JTSJ)], it 
follows that in the thermodynamic limit the excitation 
gaps (which determine the strength of the FQHE) will 
approach those of the "b" contribution alone (which is 
independent of R/Ib)- As a result, the FQHE excita- 
tion energy gaps for a graphene sample with zigzag edges 
(and in particular for v = 1/3) are 1/4 [see the prefac- 
tor of the bulk-only contribution in Eq. (|12p] of those 
calculated^ using boundless-graphene modeling. For 
boundless graphene, available in the literature EXD cal- 
culations in a spherical geometry reported two values for 
the FQHE energy gap: the first valued is ~ 0.1 e 2 /kIb 
for polarized electrons, while the second valued is 0.042 
e 2 / 'kIb when pseudoskyrmion effects are considered. 21 
From our disk-geometry EXD calculations, we find (in 
units of e 2 /kIb, and including the above 1/4 factor) an 
excitation energy gap for the "b" component of 0.0153 
(N = 6, M 1/3 = 45), 0.0160 (N = 7, M 1/3 = 63), 0.0145 
(N = 8, Mi/3 = 84), 0.0150 (N = 9, M 1/3 = 108), 
0.0149 (N = 10, Mi/3 = 135), and 0.0147 (N = 11, 
Mi/3 — 165) (Mi/ 3 gives the magic angular momenta 
corresponding to v = 1/3). These values extrapolate to 
an excitation energy gap of 0.0137 e 2 /kIb for N —> oo 
(with a standard error of ±0.0009). Our result is in good 
agreement with the experimentally measured^^ value 
(0.008 e 2 /kIb), with a possible added effect of residual 
disorder. 

To estimate the magnitude of the experimentally ob- 
servable insulating gap, \x(y = 0) = A_E(iV max + 1) — 
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FIG. 5: LLL reduced single-particle energies e as a function 
of the single-particle angular momentum m = |£|. Note the 
flat section of the curve that is followed by a rapid rise when 
the \ B bulk-type component reaches the edge of the graphene 
sample. The graphene sample boundary is taken at R = 151b- 
e is dimensionless (see text). 



AE(N ma x) [recall Eq. ([TBI) ], one needs to consider the 
maximum number, N max , of RWM electrons on the edge 
of the graphene disk [see point (I) above]. A rough es- 
timate is given by N max ~ 2ttR/{2o.Ib), where R is the 
radius of the graphene sample and a is a fitting param- 
eter. The variation of [i(v = 0) with B shown in Fig. 
[3£a) (calculated for R = 1 /im) is found to behave as B 1 , 
7 w 0.54; the increase of the gap with B is consistent 
with recent experiments on suspended graphene 3^ The 
evolution of the insulating gap with the size of the sys- 
tem is found to essentially saturate for large R [see Fig. 
Glib) calculated for B = 8 T]. 



B. Conditional probability distributions for the 
edge component 

The calculated electron density is azimuthally uni- 
form due to the conservation of the total angular mo- 
mentum M in the many-body EXD calculations. The 
spatial arrangement of the edge-type component of the 
N correlated electrons can be revealed in the intrinsic 
frame of reference via the two-body correlation function 7 
V(t, ro), referred to as a conditional probability distribu- 
tion (CPD). The CPD is defined as 

V(r,r ) cx (^iXJ^-rJ^-ro)!*^). (17) 

The CPDs are giving the probability of finding one elec- 
tron at position r assuming that another one is fixed at 
the point ro. The CPD for the edge ("e") component of 
the yrast cusp state having N = 6 and M = 39 is dis- 
played in Fig. 01 It is apparent that an RWM reflecting 
a single-ring arrangement is formed at the edge of the 
graphene sample. 
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V. SUMMARY 

We developed here a consistent picture of the many- 
body properties of interacting correlated Dirac electrons 
in graphene which is in correspondence with experi- 
mental findings pertaining to the incompressible FQHE 
states and insulating phase emerging at high magnetic 
fields. Key to the success of our theoretical model is 
the proper inclusion of the effect of the graphene-sample 
edge, treated here with the use of the zigzag bound- 
ary condition. This BC weights by a factor l/v2 the 
bulk-type wave functions in the two-component LLL 
Dirac spinor [Eq. ©]; accounting for a 50% depletion 
of the bulk electron density. The other component (also 
weighted by l/\/2) resides on the sample edge, leading 
to the emergence of an insulating phase. 

These findings point to the unique role of the graphene 
edge in modulating the strength of the interelectron in- 
teractions governing the correlated many-body states in 
graphene, thus providing an impetus for future explo- 
rations, including controlled treatments of the graphene- 
sample boundaries. 
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Appendix A: Solutions of the Dirac- Weyl equation 
in polar coordinates 

The derivation of the solutions of the Dirac- Weyl Eq. 
@ involves two subcases, i.e, for e ^ and e = 0. 

For £ ^ 0, the solutions of Eq. @ that are regular at 
the origin (x = 0) can be found through the substitution 
of the following form in the DW equation, 

X A = Ce~ x2 / 4 x^Mx), (Al) 

which will specify fi(x). 

For any I, and prior to invoking the zigzag boundary 
condition, the x A component is given by 

X A {x) = Ce-^^x^ x 




In Eq. (|A2|) . e = E /{hvp /Ib) is the reduced Dirac- 
electron energy, with Dp the Fermi velocity and Ib = 
(hc/eB) 1 / 2 the magnetic length. The general form of 
the \ B component depends on whether I is positive or 
negative, i.e., 

X B {x) = _ c m e -*v* x \i\-i M {ll ^ , (A3) 
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FIG. 6: The two LLL orbital components x B ( x ) (marked as 
B, top frame) and X A { X ) (marked as A, bottom frame) in 
Eq. ((3| for m = 3 (I = —3) and xr — 15 (a mixed bulk-edge 
DW spinor). The corresponding single-particle energy [see 
Eq. © and Fig. [5] is e = 4.383 x 10~ 22 . The vertical arrows at 
x — xr — R/Ib = 15 mark the position of the physical edge of 
the graphene sample. Inside the graphene sample (x < xr), 
the x B component can be approximated by the bulk-type 
orbital in Eq. (JSJ; outside the graphene sample (a; > xr), 
it develops an exponentially growing tail reflecting the fact 
that the full wave function in the entire range < x < oo is 
described via a confluent hypergeometric function exhibiting 
a zero at x = xr. The x A component is well represented by 
the edge orbital in Eq. ([6]). 



for I < — 1, and 



M (l + l-j,l + 2,?j), (A4) 

for I > 0. 

For e = 0, the two equations in the Dirac- Weyl system 
[Eq. ©] decouple. In this case, the regular at the origin 
(x = 0) solutions are: 



X A (x) = 

X B (x) = C^I'l-V*"/* (A5) 
for I < —1, and 

X A ( X ) = C A x l e x2 / 4 

X B (x) = 0. (A6) 

for I > 0. 

The interacting-electron EXD energies corresponding 
to these (e = 0) single-particle levels (taken as a basis) 
are higher than those resulting from the bulk-edge mixed 
states (e ^ 0); consequently they correspond to excited 
states lying above the yrast band. 

Appendix B: Two different types of edge states 

From early on, the concept of edge states has played a 
central role in the theory of the quantum Hall effect for 
nonrelativistic electrons #il In addition to these nonrela- 
tivistic edge states, the sublattice struture of graphene 
(relativistic, Dirac electrons) allows in high B for other 
unique types of edge states with no analog in nonrelativis- 
tic systems; in particular, the mixed (bulk-edge) states 
described by the spinor in Eq. ©. In this Appendix, we 
elaborate further on this point. 

To this effect, it is instructive first to investigate the be- 
havior of the LLL energies e ^ obtained as solutions of 
the transcendental equation (QJ. For a graphene sample 
with radius R — 151b, the variation of the single-particle 
energies with m = \l\ is shown in Fig. [5] We observe a flat 
region where e is vanishingly small (e ~ 0) transforming 
for higher angular momenta to a steeply rising branch. 
To gain insight into the nature of the single-particle lev- 
els in these two regions, we plot the two components of 
the Dirac-electron spinor for two characteristic angular 
momenta m. 

In Fig. IH1 we show the case of m = 3 (I = —3, with 
e = 4.383 x 10~ 22 ), which lies deep inside the flat re- 
gion in Fig. [5J Inside the graphene sample [x < xr), the 
\ B component can be approximated by the bulk-type or- 
bital in Eq. outside the graphene sample (x > Xr), 
it develops an exponentially growing tail reflecting the 
fact that the full wave function in the expanded range 
< x < oo is described via a confluent hypergeometric 
function exhibiting a zero at x — xr. The \ A component 
is well represented by the edge orbital in Eq. © , and as 
aforementioned it does not have an analog for nonrela- 
tivistic LLL electrons. 

We contrast this behavior by illustrating in Fig. [7] the 
spinor components corresponding to m = 90 (e = 0.349) 
which is characteristic of the steeply rising branch in Fig. 
[5) The x B component is similar now to a nonrelativistic 
edge state familiar from the theory of the integer QHE in 
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FIG. 7: The two LLL orbital components X B ( X ) (marked as 
B, top frame) and X A ( X ) (marked as A, bottom frame) in Eq. 
© for m = 90 (I = -90) and x R = 15 (a "double-edge" DW 
spinor, with both components representing edge- like orbitals). 
The corresponding single-particle energy [see Eq. and Fig. 
[S] is e — 0.349. The vertical arrows at x = xr — R/Ib = 15 
mark the position of the physical edge of the graphene sample. 
The \ B component is similar to the nonrelativistic edge states 
familiar from the theory of the integer QHE in semiconductor 
heterostructurespii The \ A component is similar to the edge 
orbital in Eq. (JfjJ). 

semiconductor heterostructuresjii while the \ A compo- 
nent exhibits again a behavior similar to the edge orbital 
in Eq. ©. 

From the above we conclude that the origin of the 
rising energy of the single-particle states with larger m 
(see Fig. [5j reflects the shift of the centroid of the \ B 
component toward the boundary while (simultaneously) 
satisfying the vanishing of this component on the sam- 
ple boundary, thus disturbing the orbital shape (see Fig. 
[7]). Note that the position of the centroid depends on 
to (compare Figs. H3 and [7|) and that (x 2 ) « 2(tt7 + 1) 
for yj (x 2 ) < xr. Working in the regime R/Ir » 1, 
the states in the flat region of the energy curve in Fig. [5j 
form an effective LLL manifold composed of mixed bulk- 
edge single-particle spinors [approximated by the Dirac- 
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electron spinor in Eq. ([3J]. 
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